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Abstract. In the present paper, recent experimental results on large scale coherent 
steady states observed in experimental von Karman flows are revisited from a statistical 
mechanics perspective. The latter is rooted on two levels of description. We first argue 
that the coherent steady states may be described as the equilibrium states of well- 
chosen lattice models, that can be used to define global properties of von Karman 
flows, such as their temperatures. The equilibrium description is then enlarged, in 
order to reinterpret a series of results about the stability of those steady states, their 
susceptibility to symmetry breaking, in the light of a deep analogy with the statistical 
theory of Ferromagnetism. 
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1. Introduction 

Describing the complexity of turbulent flows with tools from statistical mechanics is 
a long-standing dream of theoreticians. In 1949, flve years after the publication of its 
solution for the problem of phase transition in the 2D Ising model, Onsager published a 
notorious study of the statistical mechanics of the point vortex model [1], a special 
solution of the 2D Euler equations that allows to interpret the emergence of long 
lived coherent structures in terms of the pairing between vortices mutually interacting 
through a long range Coulombian potential. That Onsager chose the special case of 2D 
turbulence is probably not a coincidence: as soon as 1947, he was aware of the existence 
of the dissipative anomaly in 3D flows that precludes the use of classical equilibrium 
tools such as micro-canonical measures |2]. In other words, the non-vanishment of the 
energy flux at vanishing viscosity for 3D flows makes 3D turbulence an intrinsically far- 
from-equilibrium system, which cannot be properly approximated by crudely setting 
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the viscosity to zero in the Navier-Stokes equations. In a 2D flow, if one lets the 
viscosity go to zero while keeping the injection scale fixed, then the dissipative anomaly 
disappears. This justifies Onsager’s statistical mechanics approach [3]. Since Onsager, 
the description of 2D turbulence using statistical mechanics has greatly improved. 
Starting from the seminal work of Kraichnan in the 1960 ’s Hi, Robert, Miller and 
Sommeria in the 1990’s P 13 E] and subsequent work from then [9], UHl IHl [IS HSj HI] , 
the use of statistical mechanics led to a description of the coherent structures that 
seems to match the observed large scale organization in experimental and numerical 2D 
turbulence dsmsiiiTiiiH]. Yet, the extension of those ideas to 3D flows has until now 
been proven unfruitful. 

An exception may come from the special case of von Karman (VK) turbulence, a 
now classical human size experiment that reaches very large Reynolds numbers of the 
order of 10® through the stirring of a fluid in between two counter-rotating propellers. 
At this value, it is generally expected that the turbulence is fully developed with a 
wide range of interacting scales [I9] . This was indeed confirmed by previous analysis of 
the turbulence properties in the middle shear layer, which evidenced scaling properties 
and intermittency corrections in agreement with other measurements in fully developed 
turbulent flows using different geometries 1201 EH Ea Eg. Some indications exist though, 
that the number of effective degrees of freedom that describe the large scales of turbulent 
VK flows is not so large : at Reynolds number around 10® — 10®, Poincare maps of 
the torque exerted by the turbulent flow on each propeller exhibit beautiful attractors 
and limit cycles [211 ES]- Those features are usually observed in dynamical systems 
with only three or four degrees of freedom — see [26] . This suggests that the system 
could in principle be efficiently described by only a few global quantities and that a 
statistical mechanics approach could be used to identify hydrodynamical analogues 
for “temperatures” or “chemical potentials”. The present paper is precisely meant 
to support this somewhat thought-provoking idea, namely that the large scales of 3D 
VK turbulent flows can be encompassed in a broad equilibrium statistical mechanics 
framework. 

The starting point of our analysis is the observation that VK turbulence is not 
isotropic. Besides, and as far as the average flow is concerned, the swirling flows obtained 
in VK devices seem to provide an example of 3D turbulence with axial symmetry. 
As previously discussed in [271 EH], axially symmetric turbulence is an intermediate 
case between 2D and 3D turbulence, for which equilibrium theories yield non trivial 
insights [29]. VK turbulence is however not axially symmetric, and the question remains 
open whether the predictions obtained using an “axi-symmetric ansatz” are relevant 
to account for the coarse-grained properties of such flows. Preliminary comparisons 
performed at large Reynolds numbers by [201 El] suggest that the steady states of 
experimental VK flows can be described in terms of a restricted set of meta-stable 
equilibria of the 3D axially symmetric Euler equations. The goal of the present paper 
is to support further this idea and show how the light of the statistical mechanics can 
be used both beyond the scope of ideal theories and beyond the scope of strictly 2D 
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flows, in order to provide a useful framework of analysis. We will evidence a deep 
analogy between the VK steady states and lattice models of Ferro-magnetism. To 
make the analogy vivid, we stick to the simplest conceptual level compatible with a 
comprehensive description of the statistical mechanics features observed in VK flows. 
The reader interested in more technical details will be referred to the other publications. 
The present paper is organized as follows. We hrst describe the experimental set-up and 
its symmetries. We briefly recall the properties of the VK steady states, and of their 
associated bifurcations. We then summarize the outcomes of several statistical theories 
associated to the “ideal axially-symmetric fluid”. Those theories are then used beyond 
their initial scope, in order to develop an analogy between the experimental VK flow 
and lattice models of ferro-magnetism. Within this analogy, the previously observed 
VK bifurcations are shown to be reminiscent of second order mean-field transitions, and 
critical exponents are measured. We conclude by a discussion of our results. 

2. Coarse-grained description of a VK flow 

2.1. Control Parameters 

The VK experimental set-up used for the present study has been thoroughly described 
in [211 [m |32] . The fluid is confined inside a cylinder of radius R = 100 mm, and forced 
through two rotating impellers of radius Rt — see Figure All the lengths will now 
be expressed in units of the cylinder radius R. The aspect ratio of our experiment is 
dehned as the distance between the inner faces of the two opposite impellers 2H = 1.8. 
Impellers are driven by two independent motors, whose frequencies /i and /2 can either 
be set equal, in order to get an exact counter-rotating regime, or set to different values 
fi 7^ / 2 - To change the viscosity, mixtures of water or glycerol with different dilution 
rates were used. 



Figure 1. Left : Sketch of the VK2 experiment. Right: Sketch of a propeller and 
definition of the oriented angle a. 
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In this paper, three main global parameters are used to characterize VK turbulence. 
(i) The Reynolds number Re = 7 r(/i + f 2 )R^v~^ — with v the fluid kinematic viscosity 
— ranges from 10^ to 5 x 10^ so that a full range of regimes can be spanned, from a 
purely laminar to a fully turbulent one. (ii) The rotation number, 9 = (/i —/ 2 )/(/i + / 2 ), 
measures the relative influence of global rotation over a typical turbulent shear frequency. 
It can be varied from —1 to +1, hereby exploring a regime of relatively weak rotation to 
shear ratio. (Hi) Finally, the torque asymmetry 7 = {Ci — 02 )/{Ci + C 2 ) measures the 
difference between the torques Ci and C 2 applied to each of the propellers. It is crucial 
to note that in the VK2 experiment, turbulence can be either generated by maintaining 
constant the frequencies or the torques applied to each of the propeller — please, see 
[SUES] for more details. 

At a finer level of description, it has been shown that the turbulence properties 
(anisotropy, fluctuations, dissipation) are influenced by the geometry of the propellers, 
viz., their non dimensional radius Rt, the oriented angle a between the blades and the 
rotation direction (see Figure [^, the heights hh and the number n of blades In the 
present paper, we consider only propellers with hb = 0.2 and focus on changes induced 
by variations of a. Those propellers are the so-called “TM60”, “TM87” and “TM73” 
propellers, whose characteristics are summarized in Table [Tj A single propeller can be 
used to propel the fluid in two opposite directions, respectively associated to the concave 
or convex face of the blades going forward. This can be accounted by a change of sign 
of the parameter a. In the sequel, we denote (—) {resp. (-I-) ) a propeller used with the 
concave {resp. convex) face of its blades going forward. 

Table [T] summarizes the parameter space that was explored in our system. 
Schematically, the influence of the propeller geometry has been explored at Re = 10^, 
6^ = 0. The Reynolds variation has been explored at 0 = 0 using the TM60 propellers 
(±). The rotation variation has been explored at Re = 10® using TM73(±), TM87(±) 
and TM60(±). The influence of the forcing type (“constant velocity” against “constant 
torque” forcing) has been studied with the TM60(-) and TM87(-) at Re = 10®. 


Propellers 

Number of blades 

a (in degrees) 

Re 

e 

All 

8 and 16 

[-90, 90] 

10® 

0 

TM60(+) 

16 

72 

[10^, 10®] 

[-1,1] 

TM60(-) 

16 

-72 

10® 

1-1.1] 

TM87(+) 

8 

72 

10® 

[-1,1] 

TM87(-) 

8 

-72 

10® 

1-1.1] 

TM73(+) 

8 

+24 

10® 

[-1,1] 

TM73(-) 

8 

-24 

10® 

1-1.1] 


Table 1. Parameter space explored in our set-up 
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2.2. Topology of the averaged steady states 

2.2.1. A gualititative description. To analyze the topology of the averaged states, 
Stereoscopic Particle Image Velocimetry (SPIV) measurements were mostly used. The 
system provides the three components of the velocity field on a 95 x 66 points grid, 
that covers a whole meridian plane of the flow through time series of about 600 to 
5000 regularly sampled values at a lOHz frequency. We also performed a few Laser 
Doppler Velocimetry measurements providing mean velocities over a 11 x 13 points 
grid covering a half meridian plane of the flow. To deal with non-dimensional velocity 
fields, those are divided by a typical “forcing velocity” defined as Vq = 27ri?(/i -|- / 2 )/ 2 . 
We write {r,(j),z) the standard cylindrical coordinates, and denote (.)o an average 
over a time series of SPIV measurements. We also use the short-hand notation 
/x) . / = Jq -1- f{r,7i,z)) to denote the spatial average of any 

quantity / = f{r,(l),z) over the PIV window. As an example. Figure shows the {r,z) 




Figure 2. Velocity fields reconstructed from SPIV measurements at Re ~ 3 x 10® 
and 9 = 0. Top: instantaneous snapshot for TM87(-|-). Bottom : after time averaging 
over 600 snapshots, for TM87(-) (left) and TM87(-|-) (right). As the velocity field is 
projected on a meridional plane that includes the rotation axis, the left part of the 
fields here corresponds to (ur, —u^,Uz)) at ^ = tt, while the right part corresponds to 
{Ur,Uff,,Uz)) at d = 0. 


dependence of the three components of a three-dimensional velocity held {ur^u^^Uz) 
reconstructed from a PIV measurement at Reynolds number approximately equal to 
3 X 10®. Although the SPIV system does not allow us to analyze in details the azimuthal 
dependence of the velocity helds, it is clear from Figure |^that the instantaneous velocity 
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field is not axially symmetric, as it is not symmetric with respect to the transformation 
r —>■ —r. Axial symmetry exists though at the level of the time averaged velocity field 
(and more generally for any quantity derived from it). At a coarse level of description, 
the topology of the average velocity fields is simple and appears to bear some kind of 
universality. It either consists of a two-cell state [|] that is symmetric with respect to the 
equatorial axis, a two-cell state with broken symmetry or a one-cell state. 


2.2.2. A more quantitative description. As first observed in [30], the axial symmetry of 
the averaged states is not just an artifact of the averaging procedure, it can also be used 
as a natural guideline to describe the topology of the steady states. As the average flow 
inside the tank is divergence free and axially symmetric [viz. symmetric with respect 
to any azimuthal change), a Helmholtz decomposion can be used to write the averaged 
velocity field as 


(u)o 


-V(r\[') X e<A. 
r 


( 1 ) 


Independently of the underlying dynamics, the azimutal component of the vorticity 
{(jj^)o = (V X (u)o) ■ is then related to the stream function through = 

—dzz^ — drr~^drr^ = —£(d'). The knowledge of ((«</,)o, ( 1 ^ 0 ) 0 ) is then sufficient 
to reconstruct the three dimensional averaged velocity field (u)o. Using such a 
decomposition of the velocity field, Monchaux et al [201 evidenced that the axially 
symmetric averaged velocity fields observed in VK set-ups were peculiar steady solutions 
of the Euler axially symmetric equations, at least in a region far from the propellers 
and the boundaries. The Euler axially symmetric equations are derived from the 3D 
(incompressible) Euler equations by considering the dynamics of a 3D velocity field 
whose cylindrical components do not depend on the azimuthal coordinate, and depend 
on r and z only — see for example [331 EZ] and Section below. Steady states of 
the axially symmetric Euler equations are obtained whenever the toroidal field ru^, 
the poloidal field u^/r and the reduced stream function -0 = rT are related through 
relations of the kind 


rUfi) = F{%lj) and u^/r — FF'{%lj)/r‘^ = Gipp) for any function F and G. (2) 


Using TM60(±) propellers for a wide range of rotation numbers, scatter plots of both the 
toroidal field {ru^)o and the poloidal field {uj^/r)o against the reduced stream function 
p) = rT showed clear functional relationship between those quantities [SHI El]- To 
understand the general trends of the topologies, it is enough to choose linear functions 
for F and G. In this linear approximation, the VK topologies are characterized by four 
constant numbers, say A, H, G and D, defined by : [§] 

{u^)o = FT -F A/r and = {B + Dr‘^)u^ + Gr. (3) 

f Since the averaged flow is meant to be axially symmetric, the description refers to the flow observed 
in either the right or left half of the PIV window. 

§ To obtain (j^ from (|^, take F{ijj) = Bip + A and G(ip) = (C + DA/B) + Dip. 
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A, B, C, and D are computed using least square fits from the scatter plots of («</,) o 
and {Ufj,)o against d'. An example is shown in Figure obtained for a large Reynolds 
number at ( 6 ^, 7 ) = (0,0). In such a case, the £t is rather good. As already noticed 
in [3^, the £t deteriorates at lower Reynolds numbers and when the rotation number 
is too high (say \0\ > 0.3). Still, Equation ([^ provides a general framework for the 
interpretation of the data. The htting procedure was carried with various propellers at 
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Figure 3. Scatter plots obtained in TM73(+) at Re = 10® for ( 0 , 7 ) = (0,0). The 
light grey dots are the data, the opaque red dots are the fits obtained from Equation 
Left: {uOo as a function of d>. Right: {coOo as a function of 


large Reynolds number in the symmetric state (7 = 0 and 6 * = 0) with LDV data, that 
provide lower resolution representations of the mean flow. The resulting values for A, 
B, C and H as a function of the propeller’s radius and angle are provided in Figure (|^. 
Because of the measurement technique, these hts are less accurate than with the PIV 
data. The LDV-measured values of A, R, C and D should therefore be here used to 
observe trends rather than providing quantitative values. We observe that in all cases, 
both A and C are vanishing. The fact that A = 0 is compatible with being hnite at 
r = 0 , while the fact that C = 0 is a consequence of the symmetry of the basic state at 
0 = 0. The value of B depends mostly on the propeller’s angle a, being positive when 
the angle is negative and negative otherwise. The absolute value of B remains fairly 
constant in between 3 and 4, regardless of the angle — except for the a = 0 case. At 
negative a, the value of D is rather low, and close to 0. Increasing a > 0 yields a linear 
decrease for D, from 0 to —20 (at a = 72°). At this value of a though, some indications 
exist, that a second branch of solutions can be found, with D = 0. This has been 
confirmed by a study of the variations of our parameters with the Reynolds number, 
using PIV data, at 0 = 0 in the symmetric state. With these better resolved data (not 
shown), we found that the coefEcient D displays a clear bi-modal behavior, with two 
branches of solution : One extends around D = 0, and the other decreases linearly with 
log Re. A closer look at the values of D for the TM60/87(+) propellers indicates that 
the negative branch of solutions corresponds to a branch that connects continuously 
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from a two-cell to a one-cell solution. As for the other coefficients, we found that B is 
rather insensitive to the Reynolds number while both the A and C coefficients remain 
zero, at any Reynolds number, in agreement with the previously described regularity 
and symmetry properties. 




a 


Figure 4. The constant d, B, C, D as defined by Equation ([^ to characterize the VK 
topologies . Left: A (orange squares) and B (blue circles) versus the angle a; Right: C 
(orange squares) and D (blue circles) versus a. The size of the symbol is proportional 
to the impeller’s radius Rt = 0.925 (big), Rt = 0.75 (medium) or Rt = 0.5 (small). 


2.3. Transitions between the various topologies 

When the forcing is fully symmetric {9 = 0 and 7 = 0), all the impellers that we have 
tested yield a symmetric two-cell state similar to the one depicted in Figure Similarly, 
when the forcing is clearly non-symmetric (I 7 I ~ 1 or \6\ ~ 1 ), the fluid is globally in 
rotation and the average state is one-cell. Yet, the nature of the transition between the 
symmetric and the non-symmetric states does strongly depend on the geometry of the 
impellers. On the one hand, the use of low curvature impellers {a > —30°) yields a 
continuous transition that occurs via a sequence of increasingly non-symmetric two-cell 
states, with one cell becoming larger at the expense of the other. The sharpness of 
the transition can be characterized throughout the use of susceptibility coefficients (see 
Section]^, analogous to the magnetic susceptibilities in the theory of ferro-magnetism. 
Cortet et al [3ll |35] observed that those susceptibilities diverge at a hnite turbulent 
Reynolds number {Re ~ 4 x 10^), a feature clearly reminiscent of second-order phase 
transition in statistical physics. On the other hand, the use of high curvature impellers 
{a < —30°) yields an abrupt change, that gives rise to multi-stability between the two¬ 
cell symmetric state, and one of the two one-cell states (symmetric to each other with 
respect to the equatorial axis) m- As a result, a hysteresis cycle for 7 is described 
when the rotation number 6 is used as a control parameter, and cycled from 1 to — 1 
and back, over a given time scale. Increasing the curvature of the blades increase both 
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the width and the height of the hysteresis cycle that the system describes in the 0 — 7 
plane|23]. When the torqne nnmber 7 is nsed as a control parameter [i.e. when the 
system is forced at constant torqne rather than at constant velocity), the hysteresis 
cycle is regnlarized [36l |25] — see Fignre Becanse of the apparent simplicity and 



Figure 5. Variation of the topology of the stationary state, as a function of the 
forcing type. The symbols trace the torque asymmetry 7 versus the averaged rotation 
number 9. Left: for TM87(-) at constant torque (circles) and constant speed forcing 
(triangles). Right: forTM87(+) (red circles) , TM73(+) (yellow squares) and TM73(-) 
(green stars) at constant speed forcing. In insert, the corresponding topologies of the 
velocity field are shown. 


nniversality of the steady states, there is a good hope that a global nnderstanding of 
the steady states can be provided throngh general argnments based on symmetries and 
conservation laws. This is precisely the ontcome of statistical physics. In the seqnel, we 
try to explain the topology of the axially symmetric mean velocity held and explain their 
stability as a fnnction of the control parameters Re, 9, 7 , a, nsing some tools borrowed 
from statistical physics. 


3. Insights from inviscid theories 

Statistical theories of tnrbulent hows have so far only be condncted in the ideal case 
of Enler eqnations with symmetries [TJ EH HSl ESI EH EH ESI EH]- ^^is section, 
we snmmarize ideas and the ontcome of the statistical theories based on the axially 
symmetric Enler eqnations. Most of the technical details are pnshed to appendix, 
in order to focns on the predictions that those theories lead to. We then nse the 
corresponding resnlts as a gnideline to nnderstand the topologies of the how inside 
the VK set-np. 
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3.1. The axially symmetric perfect fluid. 

Inside an axially symmetric domain T>, an axially symmetric fluid is described in terms 
of a 3D velocity field u , whose components in cylindrical coordinates {ur,u^,Uz) do 
not depend on the azimuthal coordinate 0 but on r and z only. The evolution equation 
of a perfect axially symmetric fluid is therefore obtained by setting = 0 in the 
incompressible Euler equations 

dfU + u ■ Vu + Vp = 0 and V ■ u = 0. (4) 

Rather than prescribing the 3D velocity held, a convenient description of an axially 
symmetric flow can be achieved in terms of only two fields : i) the azimuthal velocity 
(later named “toroidal held”), and ii) the azimuthal vorticity held = (V x u) ■ 
(alternatively named “poloidal held”). The entire 3D velocity held can then be 
reconstructed by using the incompressibility condition (V ■ u = 0) that allows the 
following Helmholtz decomposition for u: 

u = + ^V(rT) X e 0 , (5) 

where T is the stream function, deduced from the azimuthal vorticity by the relation : 

u^ = - + dr^drU^ T = (6) 

We also define if = rT the reduced stream function. To invert the differential operator 
C, one needs to work with specified boundary conditions. We will assume here that 
is vanishing at the boundaries, a condition that comes for the impenetrability condition 
for the velocity held on the outer walls. Finally, for a cylindrical domain T> with height 
2h and radius R, we will later use the short-hand notation = (l/(hi?^)) ff^rdr dz 
later. This simple geometry will serve as a guideline for the statistical theories described 
in the remainder of the section. 


3.2. Analogy with a Ginzburg-Landau theory 


A basic input of the statistical physics of axially symmetric flows is the existence of 
conserved global quantities in the inviscid, unforced limit — at least provided that the 
perfect fluid can be considered to remain “sufficiently regular”. For instance, it is well 
known that the Euler dynamics (^ preserve both the kinetic energy E = 1/2 and 
the Helicity H = u ■ (V x u) [12]. In terms of and these quantities read : 


E 


1 

2 



ul + 


and H 





(7) 


Axially symmetric velocity fields are three-dimensionnal fields which are invariant with 
respect to the azimuthal coordinate. This spatial continuous symmetry around the axis 
makes the axially symmetric Euler dynamics further preserve two infinite families of 
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global quantities, viz. the toroidal Casimirs Cf and the generalized helicities Hg (see 
[33l [27] h This is a consequence of Noether’s theorem and of the fact that the Euler 
dynamics is Hamiltonian HDIig. The additional invariants read 


Cf = f {ruf,) and Hg = g {ru^) uj^/r for any regular enough / and g. ( 8 ) 

Jv Jv 

Those quantities allow to interprete thermodynamically the linear description of the VK 
topologies provided by Equations (§. @ and ([^. The idea is to use the invariants of 
Equation ([^ to build a “free energy” functional, that can be thought of as an analog of 
a Landau-Ginzburg functional [301 EH ES]- To stick to a linear level of description, one 
uses a quadratic form for the functions / and a linear form for g, viz., f{x) = ax^ — fix 
and g{x) = ho: — 7 , with a, fi, h, 7 yet unprescribed scalars. The free energy then reads 

U 4 ,] = / (/3/2 + ar‘^)ul + hu^oj^ + /2 - ^uj^/r - firuf,. (9) 

Jv 

If the fluid inside the VK set-up was both perfectly axially symmetric and inviscid, then 
we could expect the local minimizers of to play a peculiar role. Indeed, the reader 
familiar with dynamical systems may have already recognized that the free energy 0 
is an Arnold function relevant for axially symmetric perfect fluid, whose minima (if 
any) provide axially symmetric prohles that are formally stable — see [121133] for more 
details about the stability of inhnite dimensional dynamical systems and the stability of 
axially symmetric perfect flows in particular [j[[ The critical points of T are determined 
by the following class of axially symmetric helds (n^, : 

(/9 -I- 2ar^)M^ -1- hu^ — fir = 0 and — 7 /r = 0. (10) 


Setting A = 'f/h, B = —(5/h, C = fi/h and D = —2a/h , we exactly retrieve Equation 
(|^, which we previously used to characterize the VK topologies. This provides a clear 
connection between the steady states inside the VK tank and the steady states of the 
axially symmetric perfect fluid. However, it is easily shown that the helds (n^, 
that satisfy ( 10 ) do not in general locally minimize the free energy unless h = 0 
(non-helical case). In this case, positive values for fd and a ensure that the {u*^,uj^ys are 
indeed minimizers of B. If h is non zero, the prohles that satisfy (10) are saddle points of 
the free energy and are therefore unstable with respect to any non-trivial perturbations. 
Were we dealing with an inviscid axially symmetric huid, would we therefore conclude 
that such “meta-stable” prohles could not be observed in the long term. It is however an 
experimental fact, that those prohles (with non zero “h”) are relevant to approximate 
those observed in a VK how PI. 

To connect the VK topologies with the meta-stable prohles of Equation (10), some 
crude identihcations need to be made : the equilibrium helds (m^,^^) of the inviscid 
theory with the averaged PIV helds ( 1 ^ 0 ) 0 ), and the axially symmetric domain 


II The most general free energy for the axially symmetric dynamics is the functional E + Cf + Hg for 
prescribed / and g. Equation is a particular case. 
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V with the PIV measurement domain Vp„. The VK analogs of the axially symmetric 
Casimirs are the Angular Momentum I, the Circulation P, the Helicity H, the toroidal 
energy T, and the poloidal energy P, defined as 


I = 




{ru^)c 


P = 


'v„ 


{^4>)z 


H = 




O (^(/)) c 


( 11 ) 


T =(1/2) / (m 0 )o^ and P = (1/2) / {ur)o‘^ + {uz)c 






Equations (10) and ([^ actually provide a good approximation of the steady states of 
the VK set-up.. Figure ([^ shows the comparison between steady states obtained using 
TM73(-I-) propeller and solutions predicted by inviscid thermodynamics. To obtain it, 
we solved the equations (|^ with the SPIV-measured values of the constants A, B, C 
and D — and appropriate boundary conditions [3S]- To make the computation easier, 
we approximated the quantity B + Dr^ by a constant iPeff- We will later refer to this 
approximation as a “Beltrami” approximation. 


3.3. The statistical mechanics perspective. 

3.3.1. Analogy with a lattice model : a qualitative description. Since the topology of 
VK flows can be retrieved from a simple thermodynamic argument, one may wonder 
about the practical use of a full statistical theory. However, we also noticed that the VK 
steady prohles did not match any axially symmetric genuine equilibrium, but rather a 
class of axially symmetric “meta-stable” equilibria. One aim of statistical mechanics is 
to understand this discrepancy , and highlight important distinctions between VK flows 
and truly axially symmetric fluids. In particular, we will argue that although the VK 
prohles are in principles meta-equilibrium ones, they can still be interpreted as prohles 
that maximize a suitably-dehned conhguration entropy. 

The statistical origin of the coarse-grained steady states can be qualitatively 
intuited, by looking at the averaged signs of the azimuthal held for nearly symmetric 
forcing (7 ~ 0 and 6 ^ ~ 0). In a simplihed interpretation of the VK experiment, one may 
want to think about the signs of the instantaneous azimuthal velocity held ( as measured 
on each position of SPIV grid ) as a “spin” that could take either a -|-1 or — 1 value. In the 
sequel we drop the precautional commas “” and call it a spin. Then, each propeller could 
be thought of as a statistical reservoir of and ensuring the numbers ofand 
to remain steady. The scatter plot of the average azimuthal velocity sign against the 
stream function makes a hyperbolic tangent law emerge whatever the Reynolds number 
— see Figure At a qualitative level, the tanh law is reminiscent of the sinh laws 
observed in decaying 2D turbulence |13] , or to some particular “two-level discrete case” 


^ The equation to be solved is then simply obtained by combining the two lines of Equation (j^ — 
or equivalently (10) — and recalling that = —£4'. It reads : = KesBijj + AKes/r + Cr 

and is solved with 'I' = 0 at the boundaries ( “vanishing boundary conditions”) or by prescribing that 
u^{r,z = ±H) describe a solid rotation at frequencies fi and A (“VK boundary conditions” : see [3S] 
for more details). 
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Figure 6. Comparison between the velocity fields obtained in the VK experiment 
with TM73(+) propellers (top) and the Beltrami approximation (bottom) for “VK 
boundary conditions” at 0 = 0 (left), 6 = 0.05 (Middle) and 9 = 0.09 (right). The 
height and radius of the cylinders in the numerics are those of the real VK set-up, 
namely 2H = 1.8 and R = 1. 


found in the statistical theories developed for 2D inviscid flows in [TJ El [37] . It is also 
reminiscent of the mean-held closure equations that appear in the study of long-range 
lattice models of ferromagnetism. Consider for example the Curie-Weiss model, one 
of the simplest lattice model of ferromagnet that can be treated analytically H. and 
whose Hamiltonian reads: 

^ V N 

with Si = ±1. (12) 

i=l j^i 1=1 

In the thermodynamic limit {N ^ 1), it is well-known that the canonical free energy 
per site at temperature of the Curie Weiss model can be written as an inhmum over 
all the possible values of the magnetization /i = (1/iV) namely 

/(/3) = inf - logcosh/3(/i-b h)| . (13) 

The inhmum is reached for the value of the magnetization that satishes a self-consistent 
tanh law, namely /i = tanh/d(p + h). In the Curie-Weiss model, the tanh law emerges 
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as a consequence of the interactions being long range and the mean-field approximation 
being exact. In fact, we show in the sequel that the tanh relation between the stream 
function and the averaged azimuthal velocity spins in the VK set-up has some origin in 
this analogy with the Curie-Weiss model. 




Figure 7. The averaged azimutal velocity sign scatter plotted against the reduced 
stream function ip = rd>, at 0 = 0 ,7 = 0 and Re ~ 3 x 10'^ for TM60(-|-). ip has been 
normalized so that it ranges from —1 to +1. The light grey dots are the data, obtained 
after a time SPIV averaging. The red crosses are obtained by further averaging the 
signs on blocks of 4^ contiguous SPIV-lattice points. The black line indicates a tanh 
law. 


3.3.2. Analogy with a lattice model : the Euler perspective. The statistical mechanics 
of inviscid fluids, as developed by Robert, Sommeria and Miller for the 2D case, 
aims to determine the coarse-grained configuration (if any) that a perfect fluid is the 
most likely to adopt, were it described in terms of an equilibrium ensemble, be it micro 
canonical, canonical, grand-canonical. The coarse-grained description can usually be 
achieved in terms of a “macro-state probability field” p. In the case of axially symmetric 
fluids, the latter is defined as 

Pr(cr, 0 = Proba(rn<^ = a and ^ in the vicinity of r) (14) 

As the number of degrees of freedom of the underlying axially-symmetric Euler equations 
is formally infinite, the statistical ensembles and corresponding macro-state probability 
fields need to be defined through an appropriate discretization and coarse graining, 
whose details are here omitted — see m for the axially symmetric case and [371115] for 
the 2D case. This construction provides a useful analogy with standard lattice models of 
Ferro-magnetism. Indeed, the statistical ensemble that we wish to compute can now be 
thought of as the thermodynamic spin-wave limit of a finite-size bi-dimensional lattice 
model of interacting spins S, whose Hamiltonian is prescribed by the specific shape 
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of the inviscid dynamical invariants given by Eqnation (|^. In the axially symmetric 
case, the spins Sjj = {(Tijj^ij) are two-degrees-of-freedom objects that represent the 
micro-scale valne of the pair {rU(p,oj^/r) at position (1 <i,j < N). The analogy is 
clearer in the particnlar sitnation where the at/s can take only two valnes say a± = ±1 
(which would be the natural assumption if we wanted for example to study the 
statistics of the signs of the a^s only). In this case, the appropriate micro-canonical 
lattice model interaction Hamiltonian per spin is 

TL/N'^ = Jkiiijiki with J a discretization of (15) 

subject to the constraint that both the toroidal magnetization I = {l/N"^) 
the poloidal conditional magnetizations r± = (1/iV^) Xliyo- =±'^o remain constant. It 
however turns out that the statistical ensembles related to this “axially-symmetric lattice 
model” are ill-dehned in the thermodynamic limit [23 ES]. The problem comes from 
the poloidal degrees of freedoms, namely the ^p ’s, which are not sufficiently constrained 
by the prescription that both the energy and the conditional magnetizations remain 
hnite in the limit N ^ oo. As a consequence, an ultra-violet catastrophe occurs : 
unless the averaged poloidal held is uniformly zero, the most probable conhgurations 
are those for which most of the are inhnite and mask the coherent structures. Two 
problems are now apparent. The hrst one is theoretical. It relates to the existence 
or not of equilibrium measures for the axially symmetric Euler equations. The second 
problem is more practical : how can a ill-dehned statistical theory possibly describe the 
topology of the hows observed inside the VK set-up? As previously emphasized, the 
present paper does not deal with the axially symmetric Euler equations but with VK 
turbulence. Therefore, the hrst question is clearly far beyond our present concerns. We 
now describe a strategy to deal with the second question. 

3.3.3. A phenomenological treatment for the vorticity fluctuations. 

The main idea. Over the past decade, many statistical theories have been brought 
up to describe the statistical equilibria of the ideal axially-symmetric perfect huid 
[ISl iZlETlEEl liniEg]. Those theories yield diherent outcomes. The reason stems 
from the existence of a so-called “Ultra-Violet” (UV) catastrophe in the ideal theory 
for the axisymmetric huids. In other words : due to the special nature of the inviscid 
invariants, a full equilibrium statistical theory for the axially symmetric huid cannot 
predict both a hnite non zero energy and a non trivial large-scale how. This degeneracy 
originates from the non-zero vorticity stretching term, that allows any initial vorticity 
held to grow unbounded in the limit of a vanishing viscosity. In order to write down a 
self-consistent and non-degenerate statistical theory, one therefore needs to model the 

+ Although particular, the case is easily extended to the general case where the ct^’s span a continuous 
range of values (see [A| for more details). 
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statistics of the small scales. Different Ansatz lead to different outcomes that may or 
may not be relevant to describe VK flows. In order to select the suitable Ansatz, we first 
review the gist of the theories in the next paragraph. Their outcomes are summarized 
in Tabled We then discuss the relevance of those theories to VK flows and define two 
statistical inverse temperatures. 

The Gist of the ideal theories. As just mentioned, in common to all the works on 
axially-symmetric fluid equilibria is the use of additional phenomenological assumptions 
to control the fluctuations of the poloidal degrees of freedom. Ensemble averages 
are usually computed in terms of a most probable macro state probability field that 
maximizes a macro-state entropy, determined by an appropriate use of Laplace’s 
theorem, and the saddle-point method. The idea is therefore to restrict the set of macro 
state probability fields over which the maximization is carried on. References |13 SS] 
impose that the axially symmetric flow has a vanishing toroidal field, so that an extra- 
enstrophy constraint appears for the poloidal field; References ISHEHlEn] assume that 
the poloidal field is non-fluctuating, viz., (o;^) = at any position r ; Reference |1H] 
suggests to freeze the poloidal degrees of freedom, viz write the macro-state probability 
field as Pr{o',f) = Pr{o')pj,{^\o') and prescribe the conditional probability distribution 
Pr(^|o') ; it was also suggested to impose a cut-off on the poloidal degrees of freedom, 
either used as a physical parameter [27] or as an intermediate regularization constraint 

m- 

Table presents the different outcomes that the various ansatz lead to. Results are 
here summarized for two cases. In the “Gaussian modeling”, only the quadratic and 
linear invariants enter the theory. In the two-level modeling, it is assumed that the 
toroidal degrees of freedom can only take the two values ±1. Both cases give consistent 
results and are easily extended to the general case (see the appendix). Table shows 
that the different modelings yield different outcomes. This suggests that the poloidal 
fluctuations play a crucial role in VK mixing. 


Discussion : Temperatures and VK mixing. If one indeed takes for granted that 
the axially symmetric Euler equations can be used as a guideline to understand the 
VK topologies — which we remind is far from obvious but relies on experimental 
observations —, then those statistical theories shed a qualitative light on the mixing 
processes at stake inside the VK tank. On the one hand, the assumption that the 
poloidal field is non fluctuating and fixed in space (the “frozen-r” model) so that only 
toroidal degrees of freedom do indeed mix is not satisfying. This assumption yields only 
one of the two constitutive thermodynamic equations (10). It thus does not allow to self- 
consistently determine the averaged field without a further ansatz about the poloidal 
profile. On the other hand, if one allows for too high a level of fluctuations in the 
theory (the “micro-canonical theory” ), one obtains a closed but very restricted class of 
profiles, viz., the “(h = 0)” solutions of Equation (10). Those solutions correspond to 
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Name 

Macrostate 

Gaussian 

Two-Level 

Frozen-r 

<iK-<fo(r))Pr(o-|0 

(a) = {ti- hC{r)) /D{r) 

iO = Co(r) 

(cr) = tanh(yl -t BdO 

(C) = eo(r) 

Frozen-cr 

Pr(o-) Pr(^ 

prescribed 

(a) = (C- PA(r)^(r))/D(r) 

{0=A(r){a}+B 

(a) = tanh(Ti + BiipO)) 




Gaussian 

PrON) 

with {O < -l-C)0 

(a) = (^(r) -t B(r)^(ry)/D(r) 

(C) = (G(r)+F(r)b(r))/ZI(r) 

{a) = tanh(^i + BitpO)) 

(0 = 

1 g^QA++B+i, g_^A.+B-xl, 

~*~2ic gA+H-B+'i/'-(-e^-+s-V’ 

Microcanonical 

PrON) 

with {O —>■ -too 

(cr) = A{r) 

iO = B^PO) + C 

(cr) = Ai{r) 

iO = BMC + Cl 


Table 2. Equilibria obtained using different ansatz for the azimuthal vorticity 
Greek letters denote Lagrange multipliers, Latin letters combinations of those — that 
are different in each cell. In the Gaussian modeling only the quadratic invariants are 
taken into account, while in the two-level theory it is assumed that the toroidal field 
can only take two valnes ±1. More details can be found in the appendix. 


cases where the coarse-grained toroidal held is completely decoupled from the stream 
function. Those solutions may therefore be formally relevant for an axially-symmetric 
perfect huid, but are clearly not those observed inside the VK tank. Rather, the mixing 
inside the tank seems to he in between those two extreme cases. The “frozen-a” theory 
provides such an example of mixing. The theory describes poloidal degrees of freedom 
that cannot mix independently from the toroidal ones. Observe from Table that this 
assumption leads to a tanh relation between the coarse-grained toroidal held and the 
stream-function, as apparent in Figure 

This interpretation of statistical mechanics provides a zero order approximation to the 
actual nature of the VK stirring. The diherence between the equilibrium steady states 
(the “h = 0” case) and the non-equilibrium ones is here provided through a restriction of 
phase space. This idea is similar to the notion of restricted partition functions recently 
proposed by Herbert [19] to explore the nature of inverse cascades in helical turbulence. 
Finally note that our discussion was here restricted to predictions for the averages of 
hrst-order quantities. Obviously, a statistical theory based on an “axi-symmetric ansatz” 
fails to give an accurate and quantitative view about the huctuations measured in VK 
turbulence, as those are clearly observed not to be axially symmetric — see Figure 
1^ It is however interesting to remark that the minimal “frozen-cr” theory does not 
prescribe any shape for the poloidal fluctuations. In particular, those are not prescribed 
to be Gaussian, as one could in principle expect from a statistical mean held theory 
that uses quadratic invariants as inputs. Because of the seemingly crucial role played 
by the huctuations, it is natural to investigate in more details both the poloidal and 
the toroidal huctuations inside the VK tank. To this end, we introduce two quantities. 
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which can be thought of as two toroidal and poloidal inverse temperatures, viz. /3tor and 
/Spoi- Those are defined by 


A 


-1 

tor 



and 




(16) 


These temperatures should depend on the experimental control parameters that can 
now be seen as analogous of ’’thermostats”. In Figure we show for example the 
dependence of l//3tor as a function of the angle a (left panel) and the Reynolds number 
Re (Right panel). It can be seen that the toroidal temperature increases with decreasing 
a, and increases from zero past Re ~ 10^ (the laminar/turbulent transition). The 
temperature peaks around Re = 40000 for TM60(+) impellers, and then saturates or 
slightly decreases. Both the toroidal and the poloidal temperature will allow us to 
describe VK topologies at any control parameter from a statistical perspective, that 
goes beyond the insight of the inviscid theories. 



a logio(Re) 


Figure 8. The toroidal temperature l//3tor as a function of the control parameters 
a (left) and Re (right). The symbol color codes a (defined on the left panel). The 
size of the symbol is proportional to the impeller’s radius, Rt = 0.925 (big), Rt = 0.75 
(medium) orl Rt = 0.5 (small). 


4. Statistical mechanics beyond the inviscid case : Analogy with the 
Curie-Weiss theory of Ferro-Magnetism 

A major drawback of the theories described in the previous section is their intrinsic 
rooting on an inviscid description {i.e. force-free, zero viscosity) of the VK flows. As 
such, they do not predict anything about hnite Reynolds number effects. Besides, the 
presence of forcing and dissipation, as well as the lack of instantaneous axial-symmetry 
concur to destroy the conservation of the global invariants on which the statistical 
description relies on. However, the broad description of the VK flow topologies in terms 
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of an equilibrium statistical theory allows for a thermodynamical interpretation of the 
transitions, in the spirit of the Curie-Weiss theory of Ferro-Magnetism. Here, we shall 
not specify in as much details as in the previous section the spin lattice model that we 
consider. What we retain from the analogy is the following. The VK flow can be seen 
as a lattice model, whose spins each have two components : one linked with the toroidal 
velocity the other one linked with the toroidal vorticity u,f,. Those two-component 
spins evolve under the action of both a thermostat and a symmetry breaking external 
held, that are provided by the four control parameters Re, 7 , 9, a. The ability of the spin 
to orientate itself as a function of the forcing can be traced by the local “magnetization 
vector” (rn^(r),ci;<^(r)), whose (space-time) average reads : 


M = (/,r). 


(17) 


As illustrated in Figure different shapes for the propellers give rise to different 
behavior for M, implying different “preferred orientations” for the spins. We think 
that the propensity of each spin to deviate from this orientation can be captured by 
the behavior of the two temperatures /3poi and (3toi based on toroidal vorticity and 
velocity huctuations, see Equation (16) : The more curved the propellers, the higher 
the temperatures. In the same way, increasing the Reynolds number for a given shape 
increases the temperatures. In the spirit of statistical physics, the transition from a two¬ 
cell state towards a one-cell state can be thought of as a symmetry-breaking transition. 
It is then natural to introduce a susceptibility vector x as : 


X = iXhXr) 


aM 


89 


Re,a 


(18) 


The complete analogy is summarized in Table It allows for an interpretation of the 
salient hydrodynamical observations previously observed in the VK set-up in terms of 
phase transition and critical exponents. With this analogy in mind, a new light can 
be shed on the hnite-Reynolds number “phase transition” observed in the VK set-up 
by |Ml |35]. The main signature of the phase transition is obtained by monitoring the 
behavior of the susceptibility Xi{(^ = O); order to characterize the situation of a 
weak symmetry breaking due to the external forcing. In [HIES], the control parameter 
was taken to be the logarithm of the Reynolds number, in analogy with a dehnition 
formulated by Castaing in [50] for the temperature of a turbulent flow. Here, we show 
that the phase transition can be identihed and further characterized by the fluctuations 
of and which play the role of temperatures, as suggested by the inviscid statistical 
theory. 

The phase transition is made apparent by the study of the behavior of the 
magnetization as a function of the temperature l//3tor- This is shown in Figure[^ where 
the mean magnetizations at 6 ^ = 0 for all impellers (any a) and Reynolds numbers have 
been gathered. The data collapse nicely on a well-dehned curve. Compared to a ferro¬ 
magnetic system, the behavior is here reversed. A non-zero magnetization occurs when 
the temperature is higher than a critical temperature l//d*Q,.(M) ~ 0.044 ± 0.03. The 
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Figure 9. Magnetization M = (F, I) for various rotation numbers d’s, and different 
propellers : TM87 (+)(squares), TM87(-) (diamonds), TM73 (+)(stars) and TM73 (-) 
(triangles). The color codes 9. 


“Ferro-Magnetic” Quantity 

Hydrodynamic Analog 

Name 

Spin 

(rM 0 ,a; 0 /r) 

“Beltrami Spin” 

Magnetization 

M = (/,F) 

Angular momentum and circulation 

Thermostat 


Forcing and dissipation 

Temperature 

(1/Aor, l//9pol) 

Fluctuations 

Symmetry breaking helds 

(^,7) 

Rotation and Torque numbers 

Susceptibility 

(X/,Xr) 



Table 3. Analogy between a two components spin system and the VK flow 


data then suggest, that the magnetization grows as the square root of the distance to the 
critical temperature, viz., |M| oc VVAor - l/A*or, a behavior reminiscent of standard 
lattice models with mean-held interactions between the spins. Similarly, the behavior 
of the susceptibility y seems to exhibit a divergence around a critical temperature 
V/^*or(x) — 0.044 ± 0.03 — see Figure 11 Because of the crucial role played by the 


poloidal huctuations in the statistical theories described in Section it is tempting 
to check whether some sort of huctuation-dissipation relation holds, involving both 
the poloidal and the toroidal huctuations. Under its simplest expression, a huctuation 
relation involving the susceptibility xt can be written as 


Xr f^tor I ((^ 0)0 (^ 0 ) 0 ) f^tor/f^pol- 


(19) 




To check whether such a relation holds, we have plotted in Figure 12, xt against the 
quantity /dtor/Z^poi- Indication of a linear trend is apparent. This is compatible with 
the observation made in the previous section, that the poloidal and the toroidal degrees 
of freedom are far more correlated in a VK set-up, than a purely axially symmetric 
interpretation of the dynamics would imply. This observation is in a sense compatible 
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1 / /3tor 


Figure 10. Mean magnetization as a function of temperature. Left : I (top) and 
r (bottom) as functions of l//3tor- Right : The magnetization M = y/P + as a 
function of the toroidal temperature for different propellers (TM87(+) (circles) and (-) 
(diamond) and TM73 (+)(triangles) and (-) (stars)) . In every case, the black dotted 
line indicates a fit y = ±a-y/l//3 — 1//3*, with = 0.044 ± 0.03, and various a : 
a = 2.1 for I, a = 20 for T and a = 21 for M. The black circles are magnetization 
estimates based on the height of the hysteresis cycle [35] • 


with the “frozen-cr” scenario, which tells that the poloidal degrees of freedom are 
somehow enslaved to the toroidal ones. Note that the linear trend also exists for 
Xi- In both cases, the prefactor is small (of the order of 1/1000 to 1/20000). In 
a Beltrami approximation, the ratio /9tor//Spoi in inversely proportional to the mode 
numbers, i.e. to the number of degrees of freedom [28]. This remark would therefore 
provide an interpretation of the prefactor in the fluctuation-dissipation relation as the 
inverse number of degrees of freedom. This would give an estimate for the number of 
degrees of freedom of the order of iV ~ 10^ to 10^, much smaller than the traditional 
estimate N = (that would rather give a number of the order of 10^^ ). 


5. Conclusion 

In this work, we argued that the steady topologies observed in a turbulent VK set-up, 
as well as the transitions between those, could be interpreted in terms of a statistical 
physics modeling. This allowed us to exhibit a deep analogy between the VK large scale 
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Figure 11. The susceptibilities x as a function of the temperature. The black 
dotted line is a fit &/(l//3 — 1//?*), with b = 0.04 for xi and b = 0.2 for xr and 
1//3* = 0.044 ±0.03. 
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Figure 12. Fluctuation-dissipation relation: the susceptibilities x as a function of the 
/3tor//3poi- Left: xi- Right: xr- The dotted line are linear hts, with respective slopes 
10-3 et 2 X 10"^. 

dynamics and the standard behavior of Ferro-Magnetic material. We have shown that 
the topologies of the VK steady states conld be interpreted thermodynamically, nsing as 
a gnideline an eqnilibrinm theory that relies on the axially symmetric Enler eqnations 
and an adequate modeling of the fluctuations. At hrst sight, there was no reason to 
expect that such an equilibrium theory could be relevant for the VK cases : the lack of 
axial-symmetry for the instantaneous VK dynamics, the presence of a large-scale forcing 
which precludes any kind of separation of scale working hypothesis, the non-gaussianity 
of the VK fluctuations and the UV catastrophe associated to statistical theories based 
on the Euler axially symmetric equations are for example as many reasons which should 
imply that such an approach is bound to fail. 


log(Re) I 
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However, in the present case, the forcing and dissipation that are present in the 
experiment seem to play the role of thermostats, which prescribe strong correlations 
between the degrees of freedom present in the flow. This allows for a statistical mean- 
held interpretation of the mean how, provided one plugs into the theory additional 
ansatz non prescribed by the Euler equations. In other words, the VK steady states do 
not correspond to ideal axially symmetric equilibria. Still, to hrst order, there exists 
some good indication that those can be interpreted as maximizers over an appropriate 
restricted subset of phase space of a suitably dehned conhguration entropy. This 
experimental fact provides an intuitive ground, to think about the VK experiment in 
terms of a statistical mean held theory. We have shown that this analogy goes beyond 
the mere description of the VK topologies. We showed that the transitions between 
the various topologies could be phrased in terms of phase transition and spontaneous 
symmetry breaking, shedding an unexpected analogy between VK turbulence and the 
theory of ferro-magnetism. 

The empirical coincidence between the steady turbulent states and a class of “meta¬ 
equilibria” for the Euler equation is puzzling and opens the question of how far ’’out- 
of-equilibrium” is our system. In traditional turbulence theory, the degree of non- 
equilibriumness can be quantihed by measuring the non-dimensional energy hux. In the 
VK experiment, the latter is related to the torques exerted on both propellers. One 
then observes that this quantity is minimal when the averaged topology is a two-cell 
symmetric state, obtained by forcing at positive angle a. These states appear to be 
better htted by an equilibrium theory than the other bifurcated states, obtained at 
larger 6 or for negatively curved propeller. It is also well known that the energy flux 
measured at the large scales of a turbulent flow is very close to zero. It may therefore 
well be that the large-scale topologies that we here considered match a situation where 
the deviations from equilibrium (quantihed by the hux) are sufficiently small, so that 
some kind of perturbative theory around the equilibria states may in the end be valid. 
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A. Summary of axially symmetric statistical theories. 

In this technical appendix, we explain how to obtain Table in which the outcomes 
of the diherent statistical theories based on the axially symmetric perfect fluid are 
summarized. Our point is not to discuss those theories very thoroughly, but to show 
how various ansatz for the poloidal fluctuations lead to different relations between the 
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average toroidal field and the average poloidal field. 

Because of the specific shape of the inviscid invariants related to the ideal axially 
symmetric fluid, and in particular because of the long-range, mean field nature of 
the kinetic energy, the micro canonical averages can be computed in terms of a local 
probability field : the macro state probability field. The latter is defined though a local 
discretization of the physical domain combined with a local averaging of the velocity 
field, the details of which can be found in numerous papers — see for example | 1 S 1122 ] 
for a recent exposition. In the case of the axially-symmetric perfect fluid, it is convenient 
to define the macro state probability field p as 


Pr(o', 0 = Proba(rM 0 = crandcu^ = in the vicinity of r). (A.l) 


The probability “Proba” that appears in the latter equation is a local micro canonical 
probability, which we later write {■). In what follows, we use the short-hand notation 
L,i = //r 2 dcxd^ and still denote = {1/hR^) /_\rdrdz. 

To compute the axially symmetric equilibria, one first needs to translate the 
inviscid invariants ([^ and (|^ which by definition affect the microscopic dynamics, into 
constraints for the macro-state probability fields. For axially symmetric flows, we write 
Aa- (“toroidal areas”) , Pg- (“partial circulations”) and S the macro-state constraints 
induced by the axially-symmetric Toroidal Casimirs, Helical Casimirs and the energy 
respectively, viz., 


^<xo[p] = / / Pr(c^,0<5(c^ - C^o), P^oIp] = / / aud 

Jv JcT,i Jv J(T,i 


£[p] = 



V J a,i 




(A.2) 


Another crucial macro state functional is the “macro state entropy”, defined by 


S[p\ = - 



PMA) logPr(c^,0- 


(A.3) 


V J a,i 


In principle, for a prescribed set of macro state constraints, say a{a), 7 (cr) and E, the 
micro canonical averages () can now be computed in terms of a most probable macro 
state probability field p*. This is a consequence of an extension of Laplace’s method of 
steepest descent called the Sanov theorem, p* is then obtained by maximizing the macro 


state entropy (A.3) among the macro state probability fields that satisfy the prescribed 


set of macro-state constraints, namely : 


p* = argsup{5[p]|X[p] = a{a),T4p] = 7 (a), T[p] = E} . 


(A.4) 


In principle, the optimization problem (|A.4|) can be solved explicitly in terms of Lagrange 


multipliers, and the solution glibly written as 




and Z* = 




a 


^ exp +/(a)+^p(a) ^ , with i; = 


exp 




2A 

-13Y ~ 


(A.5) 
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The Lagrange multiplier /3, and Lagrange functions / and g are determined by the 
(functional) derivatives of the partition function : 


IV 


d log Zr 

d(5 


= E, 


IV 


6 log Zr 
<5/(fTo) 


= a(cro) and 


IV 


5 log Zj. 
Sg{(To) 


= 7 ((^o)- (A. 6 ) 


The formulation (A.4)-(A.6) is however merely formal, as the optimization problem 


(A.4) is in this case ill-defined. The trouble comes from the poloidal degrees of freedom 


being in a sense not constrained enough by the macro state constraints of Equation 
(A. 2 ). The problem is apparent in the definition of the partition function (A.5) : the 
integral there involved does not in general converge. We think that this behavior 
is an avatar of the UV catastrophe encountered in the statistical theories of ideal 
3D flows. A phenomenological taming of the problem can be achieved by further 


constraining the set of macro state fields over which the optimization problem (A.4) 
is solved. This requires the use of additional ansatz, some of which we below discuss. 
In order to carry out some explicit calculations and retrieve the equations of Table ([^, 


we will consider two simplified sets of macro-constraints (A.2). (i) In the Two-Level 
modeling, we prescribe the toroidal field to be a two-level, symmetric distribution, viz., 
Q!(cr) = p5{a — 1) - 1 - (1 — p)d{a + 1). Only five constraints then remain from the set of 
constraints (A. 2 ) : the energy, two toroidal areas A±, and two partial circulations r±. 


We write f± and g± the corresponding Lagrange multipliers that appear in Equation 


(A.5). (a) In the Gaussian modeling, we relax the constraint on the toroidal distribution, 
and only prescribe its average and variance computed over the domain. Similarly, instead 
of specifying the correlations between the poloidal and the toroidal field using the partial 
circulations, we we only prescribe a value for the helicity and the circulation. We are 
then again left with five macro state constraints. As far as the optimization problem 
is concerned, this is tantamount to prescribe a quadratic / and a linear function g in 


Equation (A.9), namely 

/(cr) = —-|- fia 


and g{a) = —ha + 7 . 


(A.7) 


A.l. Ereezing the poloidal field : The “frozen-r” scenario. 


A simple but extreme way to further constrain the optimization problem (A.4) is to 


assume that the poloidal field at position r is prescribed and non fluctuating. We 


therefore solve (A.4) among the macro states of the kind: 

Pr{(T,f) = Pr{(T\f)6{^ - ^o(r))- 


(A. 8 ) 


Whatever the field fhe macro state entropy (A.3) now reads iS[p] = 

-1Pri cr|.^o) logPr(c’'|'Co)- The extremal points of 5[p] are achieved for the now correctly 
defined macro states : 


= ^exp|-/ 3 ^ + /(cT)+^o(r)^(cx) [> - ^o(r)), 


where Zf = 


exp 


a 




(A.9) 
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An easy calculation then yields : 


(cr) = / Pr{a,^)a = — h^o)/{2K + 2l3/r‘^) in the Gaussian case, 

J <7,^ 

= tanh(A + B^o{r)) in the two-level case. 


(A. 10) 


We write D{r) = 2 k + 2/3/r^), A = (/+ — /_)/2 and B = — g-)/ 2 . In both cases, 

we also obtain {^) = as indicated in Table 


A.2. Enslaving the poloidal field : the “frozen-a” scenario 

Another way to remove the UV catastrophe is to prescribe that not only the poloidal 
averages conditioned on the toroidal areas but the whole conditional distributions p(^|cr) 


are prescribed. We then choose to solve the optimization problem (A.4) over macro 
states of the kind 


Pr(c’', 0 = Pr{o')pr{^\o'), with Pr('C|c’') being prescribed a priori. 


(A.ll) 


The macro state entropy ( |A.3 ) can now be written in terms of the (prescribed) entropy 
Srfc) of the conditional distributions as 


S[p] = -I / Pr (o-) (log Pr(cr)-Sr (cr)) with Sr(cr) = - /'pr(^ I cr) log Pr(^ I cr). (A. 12) 
Jv Ja 


Dehning .^(cr) = J^^Pri^W), we hnd that the macro-states (A.ll) that solve the 
optimization problem and extremize 5[p] are those satisfying 


Pr(^) = ^ exp {sr(cr) + ^(cr)^(a) + /(cr) - fia^/ 2 r^ - / 3 ^{a)fi/ 2 } 




with Z* = exp{sr{a)+^{a)g{a) + f{a) - fia‘^/ 2 r^ - fi^{a)fi/ 2 } 


(A.13) 


In the two-level case, we write = /^^Pr(^|cr = ±), and easily obtain 


(cr) = tanh(Ai -f- Bifi) and (/) = 


exp {sr,± + ^±g± + f±- 
J2± exp {Sr,± + ^±g± + f±- / 3 ^±'ip/ 2 } 


. (A. 14) 


with Ai = (l/2)(sr,+ ^+g+ + /+ - Sr - /-P- /_) and Bi = -/?(/+ - /-)/2. 

In the gaussian modeling, one needs to specify a shape for ^(cr) and the entropies 
Sr(cr) in order to carry the calculation further. If one assumes for simplicity /(cr) = 
A(r)cr -|- B, and Sr(cr) = const., one readily obtains 


-iA + p-hB - Afifi/2 A/ \ , D 

2{. + lil2r- + Ah) {Q=A{a) + B. 


(A. 16) 


We therefore put C = 2(7^ + p — hB), D = A{k + /d/2r^ -|- Ah) to obtain the result of 
Table [2l 









CONTENTS 


28 


A.3. A regularization through dissipation : the “Gaussian” scenario 


An alternative way to give some sense to the problem (A.4) is to restrict the 


maximization procedure to the macro state helds p’s whose local poloidal fluctuations 
are everywhere bounded, namely the ones for which a prescribed A < cx) exists such that 
r(a, .^) < A. One could rather wish to plug into the integral a coefficient so 
as to make this additional macro-state constraint the exact counterpart to the poloidal 
viscous dissipative term made explicit in Section (??), but this subtlety goes beyond the 
present point. In the gaussian modeling, the partition function logarithm of Equation 


(A.5) can be computed as 


log Zr = log TT - ^ log 5(r) + — 


h/i (7 - (3'tp/2) (k /9/2r2)(7 - fdi)/2y 
4(5(r) 


(A.16) 


where the Gaussian integration requires both o and 5(r) = (k -|- /d/2r^)z/ — h^/4 to be 
positive. We deduce 




27 (k -t- /?/ 2 r^) — hjji — {k + [5/2r'^)l3'ip 


(A.17) 


4(5(r) 4(5(r) 

Hence we retrieve the expression of Table by putting A = 4/rz/ — 2 / 17 , B = h/?, D = 86, 
C = 4:'y{K -|- /3/2r^) — 2hpi and F = —2{k + /3/2r‘^)/3. In the two-level case, the partition 
function reads 


log Zr = - log 71 -log U + log /3h/2)V(4i^)_ 

2 2 T' 


(A. 18) 


We deduce 

(cr) = tanh (Ai -|- Hi' 0 ) and 


4z/ 2z/ e/±-( 9 ±-/^^/ 2 )V 2<2 ’ 


(A.19) 


with Ai = ( 1 / 2 ) (/+ - f- + {g‘i - gX)/{2v)), Bi = -l3{g- - g+)/{Aiy). To obtain the 
expression of Table one needs to develop the squares in the exponentials in both the 
numerator and the denominator of the expression for (V), to obtain A± = f± — g\/{2v) 
and B± = g±(3. 


A.^. An ideal microcanonical scenario 

Finally, in order to define a micro canonical measure for the axially symmetric Euler 
equations, the following strategy was proposed in [29] : (i) regularize the optimization 


problem (A.4) by using a temporary cutoff, (ii) solve the cutoff dependent optimization 


problem (Hi) Relax the cutoff dependency by letting the cutoff go to 00 while using an 
appropriate scaling for the Lagrange multipliers. The result of [29] can be retrieved from 


the Gaussian scenario just described and equations (A.17) and (A.19). From Equation 
(A.17), one prescribes fd (dov, h r'U h^v , 7 = 7oi^ as z/ —)■ 0 to obtain 
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From Equation (A. 17), one prescribes 
obtain 


~ (^Qiy, g± ~ f± ~ /±,o as i/ -)■ 0, to 


{a) = tanh 




and (0 


M 1 

4^2^^ e7±.o 


(A.21) 


In both cases, the average toroidal field is constant and the average poloidal field is 
linear with ip. Hence, we retrieve the expression of Table by putting A = ^/2k, 
C = 7o/ 2 - ho/i/4K, B = Bi = -/3o^/4, Ai = tanh ((/+,o - /-,o)/2) and Ci = 

(1/2) (E±S±,oe^*«)/(E± 


References 

[1] Onsager L 1949 II Nuovo Cimento (1943-1954) 6 279-287 

[2] Eyink G and Sreenivasan K 2006 Reviews of modern physics 78 87 

[3] Duchon J and Robert R 2000 Nonlinearity 13 249 

[4] Kraichnan R H 1967 Physics of fluids 10 1417 

[5] Kraichnan R H and Montgomery D 1980 Reports on Progress in Physics 43 547-619 

[6] Miller J 1990 Physical review letters 65 2137 

[7] Robert R and Sommeria J 1991 J. Fluid Mech 229 1-3 

[8] Robert R and Sommeria J 1992 Physical review letters 69 2776-2779 

[9] Michel J and Robert R 1994 Journal of Statistical Physics 77 645-666 

[10] Jordan R and Turkington B 1997 Journal of statistical physics 87 661-695 

[11] Ellis R, Jordan R, Otto P and Turkington B 2004 Communications in mathematical physics 244 

187-208 

[12] Majda A and Wang X 2006 Nonlinear dynamics and statistical theories for basic geophysical flows 

(Cambridge University Press) 

[13] Bouchet F and Corvellec M 2010 Journal of Statistical Mechanics: Theory and Experiment 2010 

P08021 

[14] Bouchet F, Laurie J and Zaboronski O 2014 arXiv preprint arXiv:1403.0216 

[15] Turkington B, Majda A, Haven K and DiBattista M 2001 Proceedings of the National Academy of 

Sciences 98 12346-12350 

[16] Bouchet F and Sommeria J 2002 Journal of Fluid Mechanics 464 165-207 

[17] Bouchet F and Venaille A 2012 Physics reports 515 227-295 

[18] Qi W and Marston J 2013 arXiv preprint arXiv:1312.2553 

[19] Frisch U 1996 Turbulence, by Uriel Frisch, pp. 310. ISBN 0521457130. Cambridge, UK: Cambridge 

University Press, January 1996. 1 

[20] Pinton J F and Labbe R 1994 Journal de Physique II 4 1461-1468 

[21] Arneodo A, Baudet C, Belin F, Benzi R, Castaing B, Chabaud B, Chavarria R, Ciliberto S, 

Camussi R, Chilla F, Dubrulle B, Gagne Y, Hebral B, Herweijer J, Marchand M, Maurer J, 
Muzy J, Naert A, Noullez A, Peinke J, Roux F, Tabeling P, de Water W V and Willaime H 
1996 EPJB 34 411 

[22] Ravelet F, Chiffaudel A and Daviaud F 2008 Journal of Fluid Mechanics 601 339-364 

[23] Herbert E, Daviaud F, Dubrulle B, Nazarenko S and Naso A 2012 EPL (Europhysics Letters) 100 

44003 

[24] Ravelet F 2005 Bifurcations globales hydrodynamiques et magnetohydrodynamiques dans un 

ecoulement de von Kdrmdn turbulent Ph.D. thesis CEA 

[25] Saint-Michel B, Dubrulle B, Ravelet F and Daviaud F 2013 Phys. Rev. Lett. Ill 234502 

[26] Eckmann J P and Ruelle D 1985 Reviews of modern physics 57 617 

[27] Leprovost N, Dubrulle B and Chavanis P H 2006 Physical Review E 73 046308 






CONTENTS 


30 


[28] Naso A, Monchaux R, Chavanis P H and Dubrulle B 2010 Physical Review E 81 066318 

[29] Thalabard S, Dubrulle B and Bouchet F 2014 Journal of Statistical Mechanics: Theory and 

Experiment 2014 P01005 

[30] Monchaux R, Ravelet F, Dubrulle B, Chiffaudel A and Daviaud F 2006 Physical review letters 96 

124502 

[31] Monchaux R 2007 Mecanique statistique et effet dynamo dans un ecoulement de von Kdrmdn 

turbulent Ph.D. thesis Universite Paris-Diderot-Paris VII 

[32] Saint-Michel B 2013 L’experience von Kdrmdn comme paradigme de la mcanique statistique hors 

quilibre Ph.D. thesis Universite Paris-Diderot-Paris VII 

[33] Szeri A and Holmes P 1988 Philosophical Transactions of the Royal Society of London. Series A, 

Mathematical and Physical Sciences 326 327-354 

[34] Cortet P P, Chiffaudel A, Daviaud F and Dubrulle B 2010 Physical review letters 105 214501 

[35] Cortet P, Herbert E, Chiffaudel A, Daviaud F, Dubrulle B and Padilla V 2011 Journal of Statistical 

Mechanics: Theory and Experiment 2011 P07012 

[36] Saint-Michel B, Daviaud F and Dubrulle B 2014 New Journal of Physics 

[37] Miller J, Welchman P B and Cross M 1992 Physical Review A 45 2328 

[38] Welchman P B and Petrich D M 2001 Physical Review Letters 86 1761 

[39] Chavanis P H 2008 Physica D: Nonlinear Phenomena 237 1998-2002 

[40] Naso A, Thalabard S, Collette G, Chavanis P H and Dubrulle B 2010 Journal of Statistical 

Mechanics: Theory and Experiment 2010 P06019 

[41] Morrison P 1998 Reviews of Modern Physics 70 467 

[42] Holm D D, Marsden J E, Ratiu T and Weinstein A 1985 Physics Reports 123 1-116 

[43] Montgomery D, Matthaeus W, Stribling W, Martinez D and Oughton S 1992 Physics of Eluids A: 

Fluid Dynamics 4 3 

[44] Nishimori H 2001 Statistical physics of spin glasses and information processing: an introduction 

(Oxford University Press) 

[45] Potters M, Vaillant T and Bouchet F 2013 Journal of Statistical Mechanics: Theory and 

Experiment 2013 P02017 

[46] Mohseni K 2001 Physics of Fluids 13 1924 

[47] Lim C C 2003 Physics of Fluids 15 478 

[48] Thalabard S 2013 Contributions to the statistical mechanics of ideal two and a half dimensional 

flows Ph.D. thesis Universite Paris Sud-Paris XI 

[49] Herbert C 2014 Physical Review E 89 013010 

[50] Castaing B 1996 Journal de Physique II 6 105-114 



